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The scattering of temperature anisotropy quadrupole by free electrons in galaxy clusters leads to a 
now well-known polarization signal in the cosmic microwave background (CMB) fluctuations. Using 
multi-frequency polarization data, one can extract the temperature quadrupole and separate it from 
the contaminant polarization associated with the kinematic quadrupole due to transverse motion of 
clusters. At low redshifts, the temperature quadrupole contains a significant contribution from the 
integrated Sachs-Wolfe effect (ISW) associated with the growth of density fluctuations. Using po¬ 
larization data from a sample of clusters over a wide range in redshift, one can statistically establish 
the presence of the ISW effect and determine the redshift dependence of the ISW contribution to 
the rms quadrupole. Given the strong dependence of the ISW effect on the background cosmology, 
the cluster polarization can eventually be used as a probe of the dark energy. 


I. INTRODUCTION 

Linear polarization of the cosmic microwave back¬ 
ground (CMB) is generated through rescattering of the 
temperature quadrupole [jjj. At the last scattering sur¬ 
face, this leads to the primordial polarization contri¬ 
bution Q, while at lower redshifts, the temperature 
quadrupole gets scattered in regions containing free elec¬ 
trons to produce a secondary polarization signal. One 
such possibility arises from reionization which generates a 
large angular scale polarization contribution correspond¬ 
ing to the horizon at the surface of reionization ji|. At 
much smaller angular scales corresponding to few arcrnin- 
utes and below, secondary polarization is generated when 
the quadrupole is scattered by free electrons in galaxy 
clusters Q. 

There are several mechanisms to generate a polariza¬ 
tion signal in clusters. However, scattering of the pri¬ 
mordial temperature anisotropy quadrupole is expected 
to dominate the total contribution A second effect 

arises from the kinematic quadrupole due to the trans¬ 
verse motion of the cluster. This contribution, fortu¬ 
nately, has a different frequency dependence than the 
contribution associated with the scattering of the pri¬ 
mordial quadrupole. 

Following initial suggestions on the cluster polariza¬ 
tion contributions [Q, recent studies have focused on de¬ 
tailed aspects related to spectral dependences ||[|, the 
statistical signature such as the angular power spectra 
of polarization j^], and the use of cluster polarization as 
a mechanism to improve the local determination of the 
temperature quadrupole [||. 

Here, we suggest an extension to the discussion in 
Ref. ||. The basic idea is to use galaxy clusters as a 
tracer of the temperature anisotropy quadrupole and to 
statistically detect its rms value. Our discussion involves 
the measurement of quadrupole as a function of redshift 
or redshift bins based on a cluster catalog. We improve 
the previous discussion in Ref. ||] by considering the con¬ 
taminant kinematic quadrupole contribution and intro¬ 


ducing the use of multi-frequency observations. We make 
a quantitative prediction for the prospects of measur¬ 
ing the redshift evolution of the temperature quadrupole 
with future CMB polarization experiments. 

The statistical reconstruction of the rms temperature 
quadrupole, as a function of redshift, is useful for sev¬ 
eral purposes. At low redshifts, the quadrupole contains 
information related to growth of structures via the in¬ 
tegrated Sachs-Wolfe (ISW) effect 0. The presence of 
the ISW effect, however, cannot be established from the 
temperature anisotropy measurements alone, due to the 
dominant Sachs-Wolfe (SW) contribution and the large 
sample variance associated with low multipoles. 

Since the ISW contribution mainly arises from low red¬ 
shifts, one expects large scale CMB contributions to be 
correlated with maps of the large scale structure [ fio|| . 
This correlation, however, is hard to establish observa- 
tionally both due to large cosmic variance and the lack 
of adequate signal-to-noise all-sky maps of tracer fields 
& One therefore needs a probe of the ISW effect that 
is not subject to a large sample variance and we propose 
the use of cluster polarization for this purpose. Using the 
polarization signal towards a sample of clusters widely 
distributed over the whole sky, one can establish the pres¬ 
ence of the ISW effect reliably. Additionally, since the 
ISW contribution is strongly sensitive to the background 
cosmology, the redshift evolution of the rms quadrupole 
can be used as a probe of the dark energy. 

The paper is organized as follows. In § [H]. we in¬ 
troduce polarization signals towards clusters associated 
with the scattering of the primordial CMB temperature 
quadrupole and with the kinematic quadrupole generated 
by transverse cluster motions and study the use of multi¬ 
frequency polarization observations as a way to separate 
these two contributions. We will then consider the recon¬ 
struction of the rms temperature quadrupole as a func¬ 
tion of redshift using polarization observations towards a 
sample of galaxy clusters. We discuss our results in § III. 
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II. POLARIZATION TOWARDS CLUSTERS 
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CMB polarization towards clusters is generated when 
the incident radiation has a nonzero quadrupole moment. 
The two dominant origins for this quadrupole moment 
are: (a) A quadrupole from primordial CMB fluctua¬ 
tions, and (b) a quadrupole from the quadratic term in 
the Doppler shift if the scattering gas has a transverse 
velocity [§f]{|. 

Towards a sufficiently large sample of galaxy clusters, 
we can write the total rms polarization contribution as |jj 

^>Tot _ j^Prim pKin 

Prim \/6 Q rms (2) 

10 ' ' T C mb 

P Kin = -^g(x)(r)(fi), 

(1) 


where r is the optical depth to scattering in each cluster: 
r = <jt J dyn e (y). For an individual cluster, the optical 
depth can be constrained through the temperature fluc¬ 
tuation associated with the Sunyaev-Zel’dovich [Q effect. 
Since we are averaging over large samples of clusters, in 
Eq. [I], we have considered the sample averaged optical 
depth, (r). fit = Vt/c gives the transverse component 
of the cluster velocity and g(x), with x = ft^/fcsTcMB, 
is the frequency dependence of the kinematic effect (dis¬ 
cussed below). 

Note that the rms quadrupole of temperature 
anisotropies, as a function of redshift, is related to the l = 
2 contribution to the CMB power spectrum as measured 
by an observer at that redshift: Qr ms (z) = hCi{z)I. 
Large angular scale observations, such as from COBE, 
constrain the variance of the temperature quadrupole 
today to be 62(2 = 0) = (27.5 ± 2.4/rK) 2 , assuming a 
power-spectrum tilt of unity |b|]. The quadrupole is 
highly uncertain due to the fact that one is limited to 
only five independent samples. Also, these observations 
do not allow a determination of the rms quadrupole at 
redshifts other than today. However, measuring the CMB 
polarization towards galaxy clusters potentially allows to 
probe the quadrupole moment at the cluster position [||. 

We can calculate the expected redshift evolution of the 
quadrupole based on the two sources of contributions, 
one at the surface of last scattering due to SW effect 
and another at late times due to the ISW effect associ¬ 
ated with a time evolving potential 0. We write these 
two contributions to the power spectrum, projected to a 
redshift 2, as 


C l (z) = Cf w (z) + C^(z) 
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Here, r z is the radial comoving distance out to redshift 
z, r rec is the distance to last scattering (z = 1100 ), and 
G(r z ) is the growth rate of linear density fluctuations 
with 5 lin {k,r z ) = G(r z )S lin (k, 0) ||. Here, 5 lin (k,r z ) is 
the fluctuation in the linear density field and P(k, r z ), in 
Eq. |2|, is its power spectrum. 

In addition to scattering via the temperature 
anisotropy quadrupole, the kinematic quadrupole asso¬ 
ciated with the transverse motion also produces a po¬ 
larization signal towards galaxy clusters To un¬ 

derstand how this polarization contribution is generated, 
first, consider electrons moving with peculiar velocity, 
/3 = v/c, relative to the rest frame defined by the CMB. 
The CMB spectral intensity in the mean electron rest 
frame is Doppler shifted 

Iv = C e a;7(l+W - 1 ’ ^ 


where x = hv/ksTcMB and y is the cosine of the angle 
between the cluster velocity and the direction of the in¬ 
cident CMB photon. When expanded in Legendre poly¬ 
nomials, the intensity distribution is 
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and contains the necessary quadrupole under which scat¬ 
tering generates a polarization contribution with Stokes- 
parameters given by Q{y') = l/lQg{x)P 2 (l — y’) 2 = 
l/10g(x)/3 2 , U(y') = 0 where y' is now the angle be¬ 
tween cluster velocity and the line of sight. 

We can write the spectral dependence associated with 
the kinematic quadrupole, relative to CMB thermal tem¬ 
perature, as 


g{x) = | coth . (5) 

Note that the spectral dependence of the kinematic 
quadrupole contribution gives a potential method to sep¬ 
arate the two contributions. This is similar to compo¬ 
nent separation suggestions in the literature as applied 
to temperature observations, such as the separation of 
the SZ thermal effect from dominant primordial fluctua¬ 
tions Q . 


A. Frequency Separation 

Here, we will consider separation based on real space 
data instead of the multipolar analysis considered for the 
temperature ju|. We will briefly discuss the separation 
problem in the presence of noise. We refer the reader 
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to Ref. |l(| for further details. In the present case, the 
problem involves the extraction of the polarization as¬ 
sociated with the primordial quadrupole, with the kine¬ 
matic quadrupole considered as a source of confusion. 
We can write the total polarization towards a sample of 
clusters within a redshift bin as 

P = Q + K + N, (6) 

where Q is the contribution due to the primordial 
quadrupole, K is the contribution due to the kinematic 
quadrupole and N characterizes the instrumental noise. 
The two-dimensional vector field associated with polar¬ 
ization can be decomposed into two scalar quantities in¬ 
volving two orthogonal measurements, which we will sim¬ 
ply denote as ( P x ,P y ). These orthogonal quantities can 
be related to the usual Stokes parameters involving Q and 
U or the gradient and curl representation of the spin-2 
field involving E and B-modes 0- 

We assume Gaussian noise with 

(N x ) =0, 

(KN b x ) = c: b , (7) 

and similarly for the y-component. Here, a, b = 1,..., N c h 
labels the different frequency channels of the experiment. 

In Fig. [I], we illustrate the frequency dependence of 
the primordial and kinematic quadrupole contributions 
using all-sky maps of the expected polarization. In these 
plots, each polarization vector, shown as lines, should 
be considered as a representation of the polarization to¬ 
wards a cluster at that location. The polarization pat¬ 
tern due to temperature quadrupole scattering is uni¬ 
form and traces the underlying temperature quadrupole 
distribution. In describing the kinematic quadrupole we 
assume, for illustration purposes, a transverse velocity 
field with (/3 t ) ~ 10 -3 corresponding to a velocity of 300 
km sec -1 . The polarization contribution due to the kine¬ 
matic quadrupole, however, is random due to that fact 
that transverse velocities are uncorrelated; the correla¬ 
tion length of the velocity field of order 60 Mpc corre¬ 
lates velocities within regions of 1° when projected to 
a redshift of order unity. Such a correlation scale is 
much smaller than the all-sky maps we have considered 
in Fig. |. 

As shown in Fig. 1, even at high frequencies where the 
kinematic quadrupole is increased due to the spectral de¬ 
pendence, the primordial polarization still dominates the 
total contribution. This is also true when allowing for the 
fact that we may have underestimated the transverse ve¬ 
locity field; while the maximum primordial quadrupole 
related polarization remains at 4.9r /iK, the polariza¬ 
tion associated with the kinematic quadrupole scales as 
0.27y(x)(/3 t /0.001) 2 r gK. To get a comparable contribu¬ 
tion at frequencies of order few 100 GHz with g(x) ~ 
few, the transverse velocity field should involve magni¬ 
tudes of order 1000 km sec -1 ; such a high velocity is 
unlikely to be present in the large scale structure. The 


fact that primordial polarization dominates over the kine¬ 
matic contribution even at high frequencies is useful for 
the present study since the reconstruction of the primor¬ 
dial quadrupole is not likely to be confusion limited. 

To find the best-fit estimator V Q associated with Q, 
and its variance, we calculate, for say one of the compo¬ 
nents of polarization, 

M p -”-S 9 ' n )) =o - <8) 


Here, i = 0 denotes the Q contribution and i = 1 denotes 
the K contribution. g 1 defines the spectral dependence of 
both effects; by definition g° = 1 as we consider variation 
of the kinematic quadrupole contribution relative to the 
thermal CMB. g 1 is given by Eq. 0. 

Following Ref. [|u| , we find the foreground degradation 
factor (FDF) as 


FDF 


1 


1/2 
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where the dot product is defined as 


N ch 

A-B = <t$ 1 2 A * C a b l B b , (10) 

a,b— 1 


with variance in the absence of the kinematic quadrupole 
contribution given by 
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Finally, the variance in the presence of the kinematic 
quadrupole contribution is given by 

4 X = (FDF) 2 4° x )2 . (12) 


B. Primordial Quadrupole Reconstruction 

We can calculate the expected errors in the primordial 
quadrupole reconstruction by considering the enhance¬ 
ment in the instrumental noise due to the separation of 
the contaminant kinematic quadrupole contribution. Un¬ 
der the assumption that the spectral dependence of the 
kinematic quadrupole, g{x), is known, we can write the 
final variance of the determined quadrupole polarization 
as 


<Tp = 2(FDF) 2 



(13) 


Here, for simplicity, we have taken the instrumental noise 
covariance matrix to be diagonal with variance in each of 
the polarization components given by of. The factor of 2 
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accounts for the two polarization components, ( P x ,P y ), 
with equal variance. For the case of equal variance in all 
channels, a = cq, er 2 —> 2(FDF) 2 er 2 /7V c h |l6| . 

Given the uncertainty in the primordial polarization 
determination, Eq. [uj, we can calculate the expected un¬ 
certainty in the rms quadrupole by considering Eq. 1. In 
Fig. 2, we show the redshift distribution of clusters ex¬ 
pected from planned SZ surveys as calculated using the 
Press-Sclrechter 118 mass function based on analytic ap¬ 
proaches in the literature Note that the selection 
functions of clusters in SZ surveys are essentially con¬ 
stant in mass over a wide range in redshift. An experi¬ 
ment such as the planned South Pole Telescope is likely 
to detect clusters down to a mass of 10 14 Mq |2(]]. Here, 
we assume a similar survey but roughly a factor of 2.5 
increase in sky coverage with a total area of ~ 10, 000 
deg 2 . 

We bin the redshift distribution of these clusters in 
five bins between redshifts of 0 and 5. We calculate the 
optical depth to scattering in each of these clusters by 
assuming constant gas mass fractions and converting the 
gas mass in each cluster to a number density of electrons 
under spherical arguments. These optical depths range 
from 0.005 to 0.02 and typical polarization contributions, 
due to the temperature quadrupole, range from 0.02 to 
0.1 /iK. While one would require high sensitive observa¬ 
tions with an instrumental noise contribution below this 
level to detect the polarization contribution reliably to¬ 
wards an individual cluster, for the study suggested here, 
this is not necessary. Since we are attempting to make a 
statistical detection of the rms quadrupole, one can av¬ 
erage over large samples of clusters widely distributed on 
the sky and the final variance on the rms quadrupole can 
be improved, in principle, by the number of individual 
clusters. 

To consider the frequency separation of the cluster po¬ 
larization contributions, we assume an experiment with 
observational frequencies at 30, 100, 200 and 400 GHz 
with a final polarization sensitivity of 1 /iK in each of 
the channels. The frequency degradation factor (FDF) 
in this case is 2.07. For polarization sensitive Planck 
channels, at 147, 217 and 545 and published sensitivities 
ranging from 10 to 560 /iK, we determined a FDF of 4.1. 
In general a low frequency channel at the Rayleigh-Jeans 
limit is required to reduce the degradation in noise, due 
to frequency cleaning, to factors of 2 and below. 

Using the FDF factor for our four channel experiment 
and the redshift distribution of clusters, we can now de¬ 
termine how well the rms quadrupole can be established. 
We illustrate the expected rms quadrupole and errors 
around the top curve in Fig. 2. Here, we have taken a flat 
cosmological model with a matter density, in units of the 
critical density, = 0.35, baryon density Db = 0.05, 
dark-energy density Q x = 0.65 with an equation state 
w = —0.7, a Hubble constant h = 0.65, and spectral in¬ 
dex of n = 1 for primordial density perturbations. We 
have normalized the fluctuations to COBE |u]]. For il¬ 
lustration purposes, we also show the expected rms tem¬ 


perature quadrupole for w = 0 and w = — 1; the former 
corresponds to the case where dark-energy behaves as 
matter while the latter is the usual cosmological constant 
case. 

We also make use of the Planck cluster catalog as well 
as Planck polarization observations with expected sen¬ 
sitivities. The Planck SZ cluster catalog is expected to 
detect clusters down to a mass limit of 5 x 10 14 Mq 
p2[ . While the number counts per steradian is lowered 
(Fig. 2), one gains by the fact that the Planck catalog 
is all-sky. Still, due to poor sensitivity of polarization 
observations, the reconstruction is noise dominated and 
no useful information can be obtained by binning the 
catalog as a function of redshift. Thus, we combine all 
data to a single bin to estimate the rms temperature 
quadrupole. With Planck, one can determine the excess 
quadrupole, but it is unlikely that detailed information 
on its evolution can be established. We therefore suggest 
followup observations of the same cluster catalog either 
through individual targeted observations or as part of 
a post-Planck polarization mission that may attempt to 
detect the signature of inflationary gravitational waves. 

In estimating quadrupole reconstruction errors, we 
have ignored uncertainties associated with the determi¬ 
nation of optical depth in each cluster. In the case of 
Planck data, if the optical depth can be established at 
the level of few tens of percent or better, the final un¬ 
certainty in the temperature quadrupole reconstruction 
from polarization will be dominated by noise associated 
with polarization measurements and not in the determi¬ 
nation of individual cluster optical depth; Given current 
expectations based on cluster scaling relations, it is very 
likely that one cen determine the optical depth to tens of 
percent or better in upcoming data. 

For the suggested improved polarization experiment, 
the situation, however is different. For the polariza¬ 
tion noise there to dominate the final quadrupole recon¬ 
struction, optical depths to individual clusters should be 
known at the level of few percent or better. Such a de¬ 
termination will require an improved understanding of 
cluster gas physics so that one can reliably convert the 
observed SZ thermal and kinetic temperature decrements 
to measurements of optical depth. During the Planck era, 
with the availability of multifrequency data, we expect 
such a determination may be possible. Furthermore, po¬ 
tential studies with clusters, such as the one proposed 
here, motivate improvements to our understanding of 
cluster physics and we do not consider our inability to 
determine optical depths to be the limitating factor. 


III. DISCUSSION 

As illustrated in Fig. 2, a CMB temperature and polar¬ 
ization survey of 10,000 deg 2 down to a fixed mass limit of 
10 14 Mq can be used to reconstruct the rms temperature 
quadrupole and its evolution reliably. The temperature 
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information is used to detect and catalog clusters as well 
as to determine the optical depth to scattering. The same 
cluster sample is then followed up with multi-frequency 
polarization observations. In reconstructing the primor¬ 
dial temperature quadrupole, one is using these clusters 
as a tracer population in which each cluster individually 
provides information on the quadrupole. 

The first attempt to perform a study similar to this will 
come from Planck’s temperature and polarization obser¬ 
vations. While Planck will provide an adequate size cat¬ 
alog of clusters distributed over the whole sky, the poor 
sensitivity of polarization observations, however, limits 
one from extracting detailed information on the evolu¬ 
tion of the temperature quadrupole. The Planck catalog, 
however, can be used to potentially detect the excess of 
the temperature quadrupole rms at low redshifts due to 
the ISW effect. 

For a more detailed detection of the ISW effect, in¬ 
cluding its redshift evolution, one needs an order of mag¬ 
nitude higher sensitive polarization observations than 
Planck. An opportunity for such a study can be consid¬ 
ered with respect to planned SZ surveys with telescopes 
such as the South Pole Telescope (SPT). If these upcom¬ 
ing observations can be made sensitive to polarization, 
the cluster catalog detectable from such wide field SZ 
surveys can potentially be used for a more reliable recon¬ 
struction of the ISW signature. The advent of polariza¬ 
tion sensitive bolometers suggest that such possibilities 
are within reach over the next decade. 

The suggested reconstruction here is beyond what was 
considered in Ref. ||. Here, in addition to establishing 
the quadrupole, we have suggested the use of its evo¬ 
lution to directly observe the presence of the ISW ef¬ 
fect. Since large angular scale temperature observations 
or cross-correlations of the temperature with maps of the 
local universe cannot be used to establish the ISW effect 
reliably, cluster polarization observations provide an in¬ 
direct method for this purpose. 

The reconstruction, however, is subjected to several 
confusions. Here, we have considered the confusion asso¬ 
ciated with the kinematic quadrupole contribution and 
have suggested the use of multi-frequency polarization 
observations as a way to separate the two contributions. 
With adequate frequency coverage, this separation can 
be achieved with an increase in the variance of polariza¬ 
tion noise by a factor of up to ~ 2. 

In addition to the kinematic quadrupole polarization, 
galaxy clusters also produce polarization signals asso¬ 
ciated with two other mechanisms that generate a lo¬ 
cal quadrupole. These effects are due to the intrin¬ 
sic quadrupole associated with double scattering effects 
within a single cluster. After the first scattering, the 
CMB acquires an local anisotropy both due to thermal 
and kinematic effects. Second scattering then produces 
polarization. The polarization associated with the ther¬ 
mal effect is proportional to T 1 2 3 4 kT e /m e c 2 . The kinematic 
double scattering induces a polarization signal propor¬ 
tional to t 2 /3. The frequency dependence of the polariza¬ 


tion will be the same as that of the thermal or kinematic 
effect, respectively. Both second scattering effects are 
expected to be small due to their dependence on r 2 . Ad¬ 
ditionally, for double scattering with a kinematic intrinsic 
quadrupole, when averaging over large samples, one ex¬ 
pects (3 « 0 due to the randomness associated with the 
large scale velocity field. This will further suppress the 
contaminant signal. Including the double scattering asso¬ 
ciated with a thermal intrinsic quadrupole leads to biases 
in the rms quadrupole estimate of the order of a few per¬ 
cent in /iK. For unbiased estimates of the quadrupole, 
one can potentially account for this contaminant with 
prior knowledge on r and T e . 

Besides establishing the presence of the ISW effect re¬ 
liably, the reconstruction is also useful for cosmological 
purposes. As illustrated in Fig. 2, the ISW contribution, 
and thus the rms quadrupole, is strongly sensitive to ba¬ 
sic properties of the dark energy such as its equation of 
state. Though not considered in this paper, the ISW con¬ 
tribution also varies significantly with clustering aspects 
of the dark energy component beyond the usual smooth 
assumption, such as in generalized dark matter models 
of Ref. ^3| . Considering a Fisher-matrix approach to the 
measurements shown in Fig. 2, assuming a cosmological 
model with the fiducial parameters listed earlier, we de¬ 
termined that one can extract the equation of state of 
the dark energy, w with an uncertainty of ~ 0.1. Beyond 
estimates of cosmology, we suggest that in the future a 
reconstruction such as the one discussed here should be 
considered to establish the presence of the ISW effect and 
its evolution. 
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FIG. 1. The CMB polarization due to galaxy clusters. The 
polarization vectors give a representation of the expected po¬ 
larization from a cluster at the corresponding location in the 
sky. The scale is such that the maximum length of the line 
corresponds to a polarization of 4.9 r fi K. 

The background color represents the temperature 
quadrupole. The top plot shows the resulting polar¬ 
ization contribution due to scattering of the primordial 
temperature quadrupole alone. The middle map show 
the total polarization contribution (primordial and kine¬ 
matic) at 600 GHz; the total is clearly dominated by the 
primordial contribution. To highlight the kinematic po¬ 
larization, in the bottom map, we arbitrarily increase its 
amplitude by a factor of 100. Note that even at high fre¬ 
quencies where the kinematic quadrupole is increased due 
to its spectral dependence, given in Eq. |^, the primordial 
polarization still dominates the total contribution. 
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FIG. 2. Top: The expected redshift distribution of galaxy 
clusters in SZ surveys. The upcoming surveys, such as the 
one planned 4000 deg 2 observations with the South Pole Tele¬ 
scope, can probe mass limits down to 10 14 Mg and with 
Planck surveyor down to 5 x 10 14 Mg. Bottom: The rms tem¬ 
perature quadrupole as a function of redshift. The error bars 
on the top curve illustrate the expected uncertainties in the 
reconstruction of rms temperature quadrupole from followup 
multi-frequency polarization observations of galaxy clusters 
with a catalog down to 10 14 Mq and covering 10,000 deg 2 
(see text for details). The large single error on the middle 
curve illustrate the expected uncertainty with the whole-sky 
Planck cluster catalog down to a mass limit of 5 x 10 14 Mq 
and using Planck polarization observations of same clusters. 
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